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Abstract — A simple method to generate a two-dimensional 
binary grid pattern, which allows for absolute and accurate 
self-location in a finite planar region, is proposed. The pattern 
encodes position information in a local way so that reading a 
small number of its black or white pixels at any place provides 
sufficient data from which the location can be decoded both 
efficiently and robustly. 

Index Terms — de Bruijn sequences, M-sequences, self-location 
patterns 



I. Introduction 

TAKE a blindfolded man on a random one-hour walk 
around town and then remove his blindfold. How will 
he know where he is? He has several options, based on the 
information he can gather. The man could carefully count his 
steps and take note of every turn during the blindfolded walk to 
know his location relative to the beginning of his trip. Armed 
with a navigation tool such as a sextant or GPS unit, he could 
ask the stars or the GPS satellites where he is. Lastly, he 
could simply look around for a reference, such as a street 
sign, a landmark building, or even a city map with a little 
arrow saying "You are here." 

There are numerous applications where a similar problem 
is encountered. We need to somehow measure the position of 
a mobile or movable device, using some sort of sensory input. 
Wheeled vehicles can count the turns of their wheels much 
like the man counting his steps. Similarly, many devices, from 
industrial machine stages to ball-mice, employ sensors which 
are coupled with the mechanics and count small physical steps 
of a known length, in one or more dimensions. The small 
relative position differences can be accumulated to achieve 
relative self-location to a known starting point. More recent 
technologies, such as those found in modern optical mice, use 
imaging sensors instead of mechanical encoders to estimate the 
relative motion by constantly inspecting the moving texture or 
pattern of the platform beneath them. 

Sometimes the inherent accumulating error in relative self- 
location methods, or some other reasons, make them infeasible 
or unfit for certain applications, where we would want the 
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capability to obtain instant and accurate absolute self-location. 
Given several visible landmarks of known locations, a mobile 
robot could calculate its position through a triangulation Q]. 
Alternatively, cleverly designed space fiducials (e.g., J2)), 
whose appearance changes with the angle of observation, can 
also serve for self-location. 

Much like street signs for people, there are absolute self- 
location methods that provide sufficient local information to 
the device sensors, such that the absolute positioning can be 
attained. Specifically, planar patterns have been suggested, 
where a small local sample from anywhere in the pattern pro- 
vides sufficient information for decoding the absolute position. 
A naive example could consist of a floor filled with densely 
packed miniature markings, in which the exact coordinates are 
literally inscribed inside each marking. Of course, that would 
require a high sensor resolution and character recognition 
capabilities. Indeed, there are much more efficient methods, 
which do with considerably less geometric detail in the pattern. 
Some commercial products have been utilizing this approach, 
e.g., a pen with a small imaging device in its tip, writing on 
paper with a special pattern printed on it, which allows full 
tracking of the pen position at any time J3J. 

A classic method for absolute self-location in one dimension 
is the use of de Bruijn sequences |4), 0. A de Bruijn 
sequence of order n over a given alphabet of size q is a 
cyclic sequence of length q n , which has the property that each 
possible sequence of length n of the given alphabet appears in 
it as a consecutive subsequence exactly once. Thus, sampling n 
consecutive letters somewhere in the sequence is sufficient 
for perfect positioning of the sampled subsequence within the 
sequence. Several methods for the construction of de Bruijn 
sequences have been proposed, e.g. (3), Q, (SJ. There 
is also a two-dimensional generalization, i.e., it is possible 
to construct a two-dimensional cyclic arrays in which each 
rectangular sub-array of a certain size k x n appears exactly 
once in the array. These types of arrays are called perfect maps, 
e.g. ED, HO | and they can serve as the basis for absolute self- 
location on the plane. 

Of special interest and importance in communication are 
maximal-length linear shift-register sequences known also as 
M-sequences or pseudo-noise sequences ifffl . An M-sequence 
of order n is a sequence of length 2" — 1 generated by a linear 
feedback shift-register of length n. In a cyclic sequence of this 
type, each nonzero n-tuple appears exactly once as a window 
of length n in one period of the sequence exactly once. These 
sequences have many important and desired properties ifTTI . 
llT2l . A two-dimensional generalization of such sequences 
was presented in lfl2l and are called pseudo-random arrays. 
We note also that M-sequences can be used for robust one- 
dimensional location by using their error-correction properties, 
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as analyzed in lfT3l . 

In this work we propose a simple product construction to 
generate a two-dimensional binary patterns for absolute self- 
location. The paper is organized as follows. In Section [IT] we 
present the product construction based on two sequences with 
some one-dimensional window properties. A two-dimensional 
array with optimal self-location based on sensing a cross shape 



is obtained by this construction. In Section III we prove that 
the same construction can be used for reasonably effective 
error-correction of self-location with a rectangular shape. Our 
conclusions and some interesting problems for future research 



are presented in Section III 



II. The Proposed 2D Self-Location Pattern 

Our approach for building two-dimensional arrays with self- 
location properties is based on a product of two sequences, 
one of which is a de Bruijn sequence and the other being a 
sequence in which only half of the patterns appear. 

A. Half de Bruijn Sequences 

A half de Bruijn sequence of order n is a (cyclic) sequence 
of length 2" _1 which has the property that for each possible n- 
tuple X, either X or X (the bitwise complement of X), but not 
both, appear in the sequence exactly once as a subsequence. 

There are many different ways to construct half de Bruijn 
sequences. One method, in which a half de Bruijn sequence of 
length 2 ,l_1 is generated from a de Bruijn sequence of order 
n — 1 by using the inverse of the well known mapping D, called 
the D-morphism, is described in J6). Another one is based on 
M-sequences. The following results were proved in Q. 

Theorem 1: If S is an M-sequence of order n—1 then for 
each pair of n-tuples X and X either X or X appears in S, 
with the exception of the pair which consists of the all-zero 
and all-one n-tuples. 

Corollary 1: Let S be an M-sequence of order n—1 and 
let S' be the sequence obtained from S by adding another one 
to the unique run with n — 1 ones. Then S' is a half de Bruijn 
sequence. 

B. The Construction 

For two sequences T = (t 1; . . . , tx) and S = (si, . . . , sjv) 
the product 7~<8><S is a Kx N array G in which Qij, 1 < i < K, 
1 < i < N, contains the value ti ® Sj (where ® denotes 
modulo 2 addition, also known as the XOR operator). 

Take a half de Bruijn sequence T = (ti, . . . , tx) and a 
de Bruijn sequence S = (si, . . . , s^r) of orders k and n, and 
lengths K = 2 k ~ 1 and N = 2™, respectively, and let G = 
7~<8><S. Clearly, each row in G equals either 5 or £. Similarly, 
each column equals either T or T. Thus, each row and each 
column retain their window property and can serve for self- 
location in each dimension. 

Theorem 2: Each cross shaped pattern with k vertical and 
n horizontal entries appears exactly once as a pattern in the 
array G. 

Proof: Let X be a column vector of length k and Y 
be a row vector of length n. Either X or X appears in the 



sequence T. Let X the pattern which appears. Both Y and Y 
appear in the sequence S. Crosses with vertical vector X and 
horizontal vector Y appear in G only in the portions related 

dcf def - 

to Z\ = X ®Y and Z^ = X ®Y . Moreover, the crosses in Z\ 
are complements of the crosses Z^. For each cross inside Z\ 
and Zi there are two possible assignments, depending on the 
mutual entry of the vertical and horizontal component. Each 
one of these values appears in either Z\ or Z 2 . ■ 
By Theorem |2j we can use a cross sensor array to sample 
k vertical and n horizontal pixels (with one mutual pixel) in 
order to obtain self-location. 

Corollary 2: The proposed method is optimal in terms of 
the number of sampled pixels required to achieve self-location 
with a cross of vertical length k and horizontal length n. 

Corollary 3: In the array G each sampled sub-array of size 
k x n has a unique location. 

Remark 1: Similar and more sophisticated product con- 
structions to generate arrays with low redundancy and effective 
two-dimensional error-correction capabilities, were suggested 
in various papers, e.g. fl4l . |15|. 

Remark 2: In practice, the planar domain is generally not 
cyclic. In order to retain the ability to sense all 2 k ~ 1 x 2™ 
possible locations with a sensor whose footprint is k x n pixels 
array, we extend T and S by appending their first k — 1 and 
n—1 bits, respectively, to their ends. The result is now a 
(2 k ~ 1 + k - 1) x (2" + n - 1) array. 

Example 1: An example of our proposed two-dimensional 
grid pattern can be seen in Fig. [T] It was generated using a de 
Bruijn sequence of order 4 in the horizontal axis, and a half de 
Bruijn sequence of order 5 in the vertical axis, resulting in a 
cyclic array of 16 x 16 pixels. The first column and the first row 
in the figure contain the location indexes. The second column 
and the second row contain T and S, respectively. From the bit 
values inside the grid we can decode our position. An example 
of a sensor readout is marked in the table. The sensor is a 5 
by 4 cross. The vertical readout is 10010, and the horizontal 
readout is 1000 and its unique position can be easily decoded 
from T and S. 

C. Computing the Location 

The first step in our method recovers the one-dimensional 
subsequences that correspond to the location in each dimen- 
sion. Essentially, the two-dimensional problem is now reduced 
to two independent one-dimensional decoding problems. De- 
coding the location of a subsequence in a de Bruijn sequence is 
a well-known problem. Decoding of a half de Bruijn sequence 
is done similarly. 

A classic approach of creating a de Bruijn sequence S of 
order n, requires O(n) space and 0(n • 2 n ) time to generate 
the whole sequence S |5|, [8|. This involves O(n) space and 
0(n • 2") time, with n being the order of the de Bruijn 
sequence. If running time is an issue, one could create and 
store in advance a look-up table which lists the locations of all 
subsequences. This yields O(n) time complexity, but requires 
0(n • 2") space for the table. For larger n, a more flexible 
trade-off between time and space complexity was suggested 
in |[T6l . A partial look-up table of evenly spaced locations 
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Fig. 1. The 16 X 16 product array of T Cg> S. The marked cells illustrate a readout by a cross-shaped sensor. 



called milestones is created in advance. During runtime, the 
algorithm which generates the sequence is initialized with the 
query subsequence and then iterated until one of the milestones 
is encountered. For example, this can yield 0(n ■ 2?) time 
complexity and will require 0(n • 2?) space for the table. 

In either case, implementation of the self location pro- 
cess using modern computer systems is feasible, at least 
for reasonable and practical values of n, depending on the 
application. Take n = 16 for a concrete example. It allows a 
definition of 2 16 locations, e.g., a resolution of 0.1mm over 
a range of about 6.5 meters. In the first approach it would 
take, in the worst case, about 65k simple iterations (on a 
16-bit register), which can be performed reasonably quickly 
on current modest embedded processors currently clocked at 
about tens or hundreds of Megahertz. In the second approach, 
the look-up table would consume about 128k bytes (each 
entry being a two-byte word), which is, again, a quite modest 
requirement given today's memory capabilities. 

There are more efficient methods to generate de Bruijn 
sequences [17] which can be used in case of an application 
in which k and n are much larger. The problem of decoding 
perfect maps was considered for example in 1181 . A compre- 
hensive survey on this topic was given in |19|. 



Robust self-location algorithm 

The algorithm's input is a rectangle Z = {zij : 1 < i < 
k, 1 < j < n} = (X (g) Y) © £; where X is a vertical A; -tuple 
of a vertical half de Bruijn sequence T; Y is a horizontal 
n-tuple of a horizontal de Bruijn subsequence S; and £ is a 
k x n error pattern. We assume that less than ? of the bits 
in each row of £ are ones and less than | of the bits in each 
column of £ are ones. The output is the original horizontal 
and vertical subsequences X and Y, respectively. 

• Assume that the first bit of X is b. Let D be the first row 
of Z. 

• For each row A of Z 

- if more than half of the bits of A © D are zeroes 
then the corresponding bit of X is b 

- otherwise, the corresponding bit of X is b. 

• Assign or 1 to 6 to obtain X which appears in S. 

• For each column B of Z 

- if more than half of the bits of B © X are zeroes 
then the corresponding bit in Y is a zero 

- otherwise, the corresponding bit in Y is an one. 



Fig. 2. The robust self-location algorithm. 



III. Robust Self-Location 

The cross shaped sensor is rather 'spread out', so it might 
be a disadvantage in applications. In this section we show that 
this weakness becomes an advantage for robust self-location 
when the sensor is of a rectangular shape. If we use a k x n 
pixel sensor (see Corollary [5J, we can utilize the inherent 
redundancy within the kn bits to decode the location while 
overcoming a considerable number of faulty bits readings. 
This is also a very practical choice, considering that two- 
dimensional rectangular sensor grids are the most common 
variety and are the standard choice for most applications. 

We assume that less than quarter of the bits in each row 
and less than half of the bits in each column of the input 
array are in error. As it will be shown in the sequel, this 
is a fair assumption which can account for quite strong 
noise in practical terms. The algorithm for robust self-location 
presented in Fig. [2] is a simple majority decoding. 



Theorem 3: Given a grid of size 2 k ~ 1 x 2™ and a k x n 
pixel sensor, if less than quarter of the bits in each row and 
less than half of the bits in each column are in error, then the 
algorithm accurately decodes the sensor location. 

Proof: Since the number of errors in a row is less than ~ 
it follows that two rows which were originally the same will 
agree in more than half of their bits and two complement rows 
will disagree in more than half of their bits. Therefore, the 
related bits of X will be the same or different, respectively. 
Having all the k bits of X in terms of the variable b, there 
is only one assignment of a legal fc-tuple since the vertical 
sequence is a half de Bruijn sequence. 

Having the correct vertical subsequence X and since the 
number of errors in a column is less than | it follows that 
if X agree in more than ^ bits with a column then the 
corresponding bit of Y is a zero; and if it disagree in more 
than § bits with a column then the corresponding bit of Y is 
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a one. m 

Remark 3: Decoding can be done also if more than quarter 
of the bits in some rows are in error. A slightly better condition 
would be to require that the number of distinct positions in 
error in any two rows is less than ?|. This requirement can be 
further improved. 

Similar algorithm will also work if we will exchange 
between rows and columns, or equivalently if we will consider 
a transposed array. Therefore, we can exchange our assumption 
on the number of wrong bits in a row or a column. But, 
having for example at least half of the bits wrong in a given 
column (or a given row) will cause a wrong identification of 
the original subsequences. 

Lemma 4: Given a grid of size 2 k ~ 1 x 2™ and a k x n pixel 
sensor, if at least half of the bits in one of the columns of 
a pixel sensor are in error, then we cannot ensure accurate 
decoding of the original subsequences. 

Proof: Let X and Y two n-tuples which differ only in 
the first bit. Both X and Y appears as a window of length n in 
the de Bruijn sequence S of order n. Let Z be a fc-tuple which 
appears as a window in the sequence T. The products Z ®X 
and Z®Y appears as k x n windows in the array T®S. Both 
k x n windows differ only in the first column and it would be 
impossible to distinguish between the two windows if half of 
the bits in the first column are in error. If more than half of 
the bits in the first column are in error then a wrong decoding 
of the sensor location will be made. The same arguments can 
be applied to any other column. ■ 
We note that by Lemma |4] we cannot correct |~|] or more 
random errors in a k x n array. The reason is that the array is 
highly redundant. This is quite weak from an error-correction 
point of view. But, by Theorem [3] we are able to correct 
about errors in an k x n array if less than ? errors 
occur in a row and less than | errors occur in a column. 
The reason is that redundant rows and columns are used for 
the majority decoding. This result is quite strong from error- 
correction point of view. Thus, the weakness for one type of 
errors becomes an advantage for another type of errors. 

Example 2: The 7x9 sub-array of Fig. |3]has no more than 
two errors in a row and three errors in a column. The first row 
has more than half bits in common with the 5th and the 7th 
rows. Thus the vertical pattern is bbbbbbb. Suppose that 6 = 1, 
i.e. the vertical pattern is 1000101. We now compare all of the 
columns with 1000101. If more than half of the corresponding 
bits agree, the bit in the horizontal sequence is one; otherwise 
it is a zero. Thus, the sequence is 110111001. The sub-array 
with no errors is presented in Fig. [4] 

Now, we analyze the error rates in individual bits that allow 
us to determine the probability that the position is determined 
correctly. Given a k x n rectangle in which the probability of 
each bit being correct is p independent of the other bits, we can 
determine the probability that each row satisfies the condition 
above, that less than quarter of the bits are in error. Then the 
probability that each row is satisfactory is the individual row 
probability raised to the kfi 1 power, the number of rows. For 
simplicity we assume now that k = n. 

To find the probability that all of the rows satisfy the row 
condition, we raise the probabilities P(n;p) to the power of 



Fig. 3. A 7 X 9 sub-array with errors. 



Fig. 4. The corrected sub-array. 



the number of rows. These are given in Table [I] 
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TABLE I 

Probabilities that each of the rows of an n x n window has 

FEWER THAN QUARTER OF ITS BITS IN ERROR WHEN THE PROBABILITY 
THAT EACH BIT IS CORRECT IS p. 



In order to have a probability of at least 0.99 that the row 
condition is satisfied, we need p > 0.994 for n = 8, p > 0.98 
for n = 16, p > 0.95 for n = 32, and p > 0.91 for n = 64. In 
order for the row condition to be satisfied with probability at 
least 0.999, it is sufficient that p > 0.99 for n = 16, p > 0.96 
for n = 32, and p > 0.93 for n = 64. Also, if p > 0.98 for 
n = 32 or p > 0.94 for n = 64, the row condition is satisfied 
with probability greater than 0.9999 . 

The probability that the column condition (that less than half 
the bits are in error) is not satisfied when the row condition is 
satisfied is negligible. For example, if we let Q(n;p) represent 
the probability that the column condition is not satisfied, i.e. 



Q(n;p) 



Kil 



p\\- P r 



(i) 



then we have results such as Q(16; 0.99) = 1.2 x 10~ 12 . In a 
square array, the probability that the column condition is not 

satisfied for at least one of the columns is then 1 — (1 — 1.2 x 

10-12)16 = i g x 10 -U 



IV. Conclusion and Future Research 

Implementing absolute self-location in a planar region using 
special patterns is a viable and proven approach and can solve 
a variety of technological problems. In this paper we proposed 
a solution based on robust two-dimensional arrays with a two- 
dimensional window property. The method also has a rather 
strong error-correction capability. It enables to correct errors 
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if less than quarter of the bits in a row and less than half of 
the bits in a column are in error. 

In some applications, the alignment of the sensor array to 
the grid pattern is not guaranteed. The sensor may be arbitrar- 
ily translated and rotated, so that retrieving the local bit matrix 
is not trivial. Position location of one-dimensional sequences, 
when the orientation of the subsequence is not known was 
considered in [20 1 . The solution in two-dimensional arrays is 
to sample the region at a somewhat higher resolution than 
k by n, and analyze the image in order to first estimate the 
pose of the pattern of rows and columns. Since the proposed 
pattern has a very pronounced structure consisting of identical 
or inverted rows (as well as columns), this can greatly aid in 
the task. Using an M-sequence and its complement as vertical 
and horizontal sequences in our construction can also help in 
solving of the orientation problem. A complete analysis of 
these issues is a problem for future research. 

There are many other future research problems in this area. 
Some related to our specific construction and some are to new 
possible construction methods. 

1) As indicated in Remark [3] the claim in Theorem [3] can 
be strengthened. What is the strongest claim on the error 
capability of our scheme? Do the de Bruijn sequence and 
the half de Bruijn sequence that we selected have any 
influence on this claim? 

2) How can we improve the error-correction capabilities 
of our scheme if the de Bruijn sequence and the half 
de Bruijn sequence are derived from M-sequences with 
error-correction capabilities as indicated in fl3l . 

3) The array obtained by our method can correct a limited 
number of random errors, even so we proved that the 
probability for such errors which the method cannot 
correct is negligible. Generating arrays with window 
properties which can correct large number of random 
errors is an important topic for future research. 

4) Finally, we note that a folding method for generating 
pseudo-random arrays from M-sequences was suggested 
in |[T2l . This method was subsequently generalized 
in ll2TI . Can this method be adapted also to generate 
better pseudo-random arrays which can correct random 
errors? Using the M-sequences as suggested by lfl3l for 
this purpose could be the first step in attempting to find 
an answer to such questions.. 
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